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Abstract 

In the framework of Brans-Dicke (BD) theory, the present study determines the time 
dependence of BD parameter, energy density and equation of state (EoS) parameter of the 
cosmic fluid in a universe expanding with acceleration, preceded by a phase of deceleration. 
For this purpose, a scale factor has been so chosen for the present model that the deceleration 
parameter, obtained from it, shows a signature flip with time. Considering the dark energy to 
be responsible for the entire pressure, the time evolution of energy parameters for matter and 
dark energy and the EoS parameter for dark energy have been determined. A term, 
representing interaction between matter and dark energy, has been calculated. Its negative 
value at the present time indicates conversion of matter into dark energy. It is evident from 
the present study that the nature of dependence of the scalar field upon the scale factor plays 
a very important role in governing the time evolution of the cosmological quantities studied 
here. This model has an inherent simplicity in the sense that it allows one to determine the 
time evolution of dark energy for a homogeneous and isotropic universe, without involving 
any self-interaction potential or cosmological constant in the formulation. 

Keywords: Brans-Dicke Theory, Dark Energy, EoS Parameter, Density Parameters, 
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1. Introduction 

The accelerated expansion of the universe is one of the most interesting and important 
phenomena in the field of cosmology that have been obtained through astrophysical 
observations in recent times [1-5]. An exotic form of energy, with a negative pressure, has 
been found to be responsible for the accelerated expansion of the universe. This energy is 
known as dark energy (DE). In the fields of physics and astronomy, an extensive research is 
now taking place, throughout the world, on the nature and dynamics of DE. A number of 
models have been proposed to explain this accelerated expansion of the universe, following a 
phase of deceleration. In most of the models, DE is represented by cosmological constant [6]. 

The present article is based on Brans-Dicke (BD) theory of gravitation. The BD 
theory is characterized by a scalar field ( <p ) and a dimensionless coupling parameter (co) that 
govern the dynamics of space-time geometry. It can be regarded as a natural extension of the 
general theory of relativity which is obtained in the limit of an infinite co and a constant value 



of 0 [7]. The BD theory of gravity can be regarded as one of the most important theories, 
among all prevalent alternative theories of gravitation, which have very successfully 
explained the early and late time behaviours of the universe and solved the problems of 
inflation [8]. As an extension of the original BD theory, a generalized version was proposed, 
where m is regarded as a function of the scalar field 0 [9-11]. Several models regarding the 
expanding universe have been formulated on the basis of this theory [7, 12-15]. 

In the present study, BD field equations (for a spatially flat, homogeneous and 
isotropic space-time) have been used to determine time dependence of the coupling parameter 
(to), equation of state (EoS) parameter (y) and energy density (p). Time dependence of to was 
previously studied by many groups, using scale factors that have power law dependence upon 
time [7, 14, 15]. These scale factors lead to constant deceleration parameters that are not 
likely to describe the cosmic expansion properly. According to astrophysical observations, 
the universe has made a transition from a phase of decelerated expansion to a phase of 
accelerated expansion, implying that the deceleration parameter has a dependence upon time 
[12, 16]. Taking this fact into consideration, we have used a scale factor (a) that leads to a 
time dependent deceleration parameter which changes sign from positive to its present 
negative value. For simplicity, we have assumed power law dependence of to upon 0 and 
also 0 upon a. Time dependence of p and the EoS parameter (y), in the present study, are 
consistent with the results obtained from some recent studies based on general relativity with 
anisotropic space-time [17, 18]. A theoretical model has been proposed to determine the time 
variation of density parameters (CL) for matter and dark energy. The results from this model 
are in agreement with those obtained from studies based on completely different premises 
[19, 20]. An effective interaction term (Qe//) between matter and dark energy has been 
calculated here and its time dependence has been studied. Its time variation is found to be 
consistent with studies based on some other models [21, 22]. 


2. Metric and Field Equations 


The action, in the generalized Brans-Dicke theory, is expressed as, 

+ ^0# tv d M 0d v 0 + £ m ) (01) 


In equation (1), g is the determinant of the metric tensor g^ v , L m is the Lagrangian for matter, 
R denotes the Ricci scalar, cp is the Brans-Dicke scalar field and w is a dimensionless 
coupling parameter. This parameter is regarded as a function of cp in generalized BD theory. 
The variation of action in equation (1) leads to the following field equations. 
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R^ v is the Ricci tensor and 7) iV is the energy-momentum tensor. In the above equations, a 
semicolon stands for a covariant derivative and a comma indicates an ordinary derivative 
with respect to . 

The energy-momentum tensor (T^ v ) for the cosmic constituents is given by, 
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Here, p denotes the energy density, p is the isotropic pressure, u v is the four-velocity vector 
and T is the trace of Tf. This form of the energy-momentum tensor is based on an assumption 
that the total matter-energy content of the expanding universe is a perfect fluid. 

In a co-moving coordinate system, u v — (0,0,0,1), having the characteristic of g^ v u^u v — 1. 
The line element for a homogeneous and isotropic universe, in Friedmann-Robertson-Walker 
(FRW) cosmology, can be expressed as, 

ds 2 — —dt 2 + a 2 (t) + r 2 dd 2 + r 2 sin 2 9d ^ 2 j (05) 

Where a(t) denotes the scale factor, t is the cosmic time, k is the spatial curvature 
parameter. Three spherical polar coordinates are denoted by r, 6 and / respectively. 

The following three equations are the field equations of BD theory for a universe filled with a 
perfect fluid, in the geometry of FRW space-time expressed by equation (5). These equations 
have been obtained from the field equations (2) and (3) in which 7/ v from equation (4) has 
been incorporated and the components of the metric tensor g fiv have been taken from the line 
element expressed by equation (5). 
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In these field equations, y(= P/p) is the equation of state (EoS) parameter for the cosmic 
fluid, which has been treated as a function of time in the present study. 

3. Theoretical Model 
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Combining equations (6), (7) and (8), for zero spatial curvature (k — 0), one obtains, 

d) + f2| +6--|)<u-6f4 + -)|= 0 (09) 
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In the present study we have assumed the following ansatzes for <p and a). These parameters 
have been assumed to have power-law relations with a and (p respectively. 
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The expression of <p in equation (10) has been taken from some earlier studies in this regard 
[12, 16] .The reason for choosing the above empirical relation for oo is the fact that this is 
regarded as a function of the scalar field in the generalized Brans-Dicke theory [12]. 

Using the expressions (10) and (11) in equation (9), one obtains, 
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Here, q (= —da/a 2 ) is the deceleration parameter. 

Writing m = o > 0 and q — q 0 (i.e., their values at t — t 0 ) in equation (12), we get, 
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Combining equation ( 6 ) with (10) and taking k — 0 (for flat space), one gets, 
p — (pH 2 ^3 + 3n — 

Here, H (= a/a) is the Hubble parameter. 

Replacing all parameters in equation (14) by their values at t = t 0 , one obtains, 
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Using equation (15) in equation (13) one gets, 
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Equations (15) and (16) show that both <u 0 and m depends upon the parameter n, which 
determines the time dependence of 0 , as per equation ( 10 ). 

Using equations (10) and (14) in (7), for k = 0 (zero spatial curvature) and writing ^ = 
—qH 2 one obtains the following expression of the equation of state parameter, 
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The value of the equation of state parameter at the present time is therefore given by, 
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To use the expression of co in equation (11), the values of o> () and m should be taken from 
equations (15) and (16) respectively. The value of m, required for the expressions of p and 7 , 
in equations (14) and (17) respectively, should then be taken from equation (11). The value of 
0, required for equation (14), should be taken from equation (10). 

To determine the time dependence of several cosmological parameters (0, 00 , p, y), following 
empirical expression of the scale factor has been used in the present model. 

a — a 0 £'xp[a{(t/t 0 ) /? — l)] (19) 

This scale factor has been so chosen that it generates a deceleration parameter (q — —da/a 2 ) 
which changes sign from positive to negative, as a function of time. This signature flip 
indicates a transition of the universe from a phase of decelerated expansion to a phase of 
accelerated expansion, in accordance with several recent studies based on astrophysical 
observations [12].We have taken a 0 = 1 for all calculations. 

Here a, (3 should have the same sign to ensure an increase of the scale factor with time. Using 
equation (19), the Hubble parameter ( H ) and the deceleration parameter (q) are obtained as, 
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For 0 < /? < 1 and a > 0 one finds that, q -* +oo as t -* 0 and q -* — 1 as t -> qo, showing 
clearly a signature flip of q with time. 

Taking H — H 0 and q = q 0 , at t — t 0 , the values of the constants a and /? are obtained as, 
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It is often necessary to find the evolution of a cosmological quantity as a function of the 
redshift parameter (z) where z = ^ — 1. Using equation (19) one gets the following 
expression as a relation between redshift parameter and time. 

z = 7 - 1 = Fxp[-a{(t/t 0 / - 1}] - 1 (24) 

The values of different cosmological parameters used in this article are: 
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4. Determination of Density Parameters 

The total pressure (P) of the entire matter-energy content of the universe is contributed by 
dark energy because, the whole matter content (dark matter + baryonic matter) is regarded as 
pressureless dust [12, 23]. Thus we can write, 

P = YP = Yd Pd (25) 

Here, y D is the EoS parameter for dark energy and p D is the density of dark energy. 

We propose the following empirical relationship between y D and y. 

Yd = A yt* (26) 


Here A and p are constants. 

Using equations (25) and (26) we can write, 

Q D = — = — = A-H-v (27) 

P Yd 

Here, Q D denotes the density parameter for dark energy. 

Taking Cl D — Q D0 at t = t 0 , in equation (27), we get, 
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Using equation (28) in (27) we get, 
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Using equation (29), the density parameter for matter can be expressed as, 

= p -f = ™ = 1 - = 1 - n D0 (£)~ 11 (30) 

In deriving equation (30), we have used the fact that the universe is composed mainly of 
matter and dark energy, considering all other forms of energy to be negligibly small [24, 25]. 
Using equations (29) and (30), the expressions for densities of dark energy and matter can be 
respectively written as, 
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Using equation (28) in (26), the equation-of-state parameter for dark energy (y D ) becomes, 
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Using equation (24), the expressions of Q D and Q m can be written in terms of redshift (z) as, 



The values of Q D0 is close to 0.7 according to several astrophysical observations [24, 25]. 

As per the history of evolution of density parameters of the universe, there was a time in the 
recent past when Q D — Q m = 0.5 and, the corresponding z value was lying somewhere in the 
range of 0 < z < 1, as obtained from some recent studies [19, 20]. That phase of the universe 
might be same as (or close to) the one when there was a transition from decelerated 
expansion to an accelerated expansion of the universe because it is the dark energy that is 
known to be responsible for the accelerated expansion. According to a recent study, the 
transition of the universe from a phase of decelerated expansion to its present state of 
accelerated expansion took place in the past at z = 0.6818 which was around 7.2371 x 10 9 
years ago [25]. Taking z c to denote the value of the redshift parameter at which the universe 
had Q d — Q m , we get the following expression of /r from equations (34) and (35). 
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Using equations (10), (19), (29) and (30), Q D and Q m can be expressed as functions of the 
scalar field ( <p ) by the following two equations. 
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5. Calculation of an Effective Interaction Term 

The energy conservation equation is expressed as, 

p + 3 H{p + P) = 0 (39) 

Taking the pressure P — yp — YdPd (eqn. 25) and p — p m + p D in equation (39) one gets, 

Pm + Pd + [p m + p D (1 + Yd) 1 = 0 (40) 

If it is assumed that the two entities, matter and dark energy the universe, have been 
interacting with each other, causing the generation of one of them at the cost of the other, one 
may define a parameter representing their interaction, on the basis of equation (40). The 
interaction term ( Q ) can be represented by the following equation [19, 21, 23]. 

Pm + 3 Hp m = Q (41) 

p D + 3Hp D (l + Y D ) = -Q (42) 

A negative value of Q represents a transfer of energy from the matter field to the field of dark 
energy and a positive value of Q represents conversion of dark energy into matter. It is 
evident from equations (41) and (42) that Q has a dependence upon time and, its values, 
obtained from these two equations, would not be the same. Let us denote the values of Q, 
obtained from these two equations, by Q 1 and Q 2 respectively. 

Using the relation p — p m + p D and equation (31) in (41) we get, 

Qi = (1 - n D )p - P&D + 3Hp(l - n D ) (43) 

Using equations (31) and (33) in equation (42) we get, 

Qz = -&dP ~ P&d ~ 3 Hp(Q D + y) (44) 

To estimate the difference (A Q) between Qpind Q z , we write the following expression. 

AQ = Q 1 -Q 2 =P + 3Hp(l + Y ) (45) 


To have an estimate of an effective interaction between matter and dark energy, we define a 
parameter Q e ff as the average of Qpind Q 2 . Thus we have, 
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To judge the plausibility of using Q e ff as a measure of interaction between matter and dark 
energy, one may calculate A Q/Q e ff, from the following expression (using eqns. 45, 46). 


AQ _ p + 3tfp(l+y) 

Qeff Q—£tD)p~pAD+|tfp(l-2fi D -y) 


(47) 


To determine the time dependence of the interaction term ( Q ), we have calculated p and Q D , 
based on the expressions of p and Q. D in equations (14) and (29) respectively. 
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Equations (10) and (11) have been used, along with the definitions of H and q, to calculate p 
(eqn. 48). Equation (20) has been used to calculate Cl D (eqn. 49). 

Using equations (48) and (49), in (43), (44), (46) and (47) one obtains, 
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If the value of is found to be sufficiently small, it would be quite reasonable to believe 

Qeff 

that Q e ff, which is actually an average of Qptnd Q 2 , will truly represent the interaction 
between matter and dark energy. 

In the table below, i.e. Table-1, we have listed the values of (Q e ff) and f-^-| for 

\QeffJ t=to 

different values of the parameter n. These are the values of these two quantities at the present 
time. Table-1 shows that the present value of the interaction term between matter and dark 
energy is of the order of 10 -44 and it is negative. The present value of A Q/Q e ff is so small 
(~10“ 15 ) that Q e ff, as defined by eqn. 46, is likely to describe the interaction with sufficient 
accuracy. Its negative value indicates decay of matter into dark energy at the present time. 
These values show that as n becomes more negative, (<3 e //) _ becomes more negative. It 

t £q 


clearly implies that a greater rate of change of the scalar field (<p) or gravitational constant 
(i G = 1/0) with time has certainly some connection to a larger interaction between mater and 
dark energy of the universe, causing a larger rate of generation of dark energy at the expense 
of matter. 



TABLE - 1 

Values of Q e ff and A Q/Q e ff at the present time for several values of the parameter n which 
determines the time variation of the scalar field (0). 

n 

( Qe ff\=t 0 


f A<n 

\Qeff)t = t n 

-1.60 

-1.153 x 10“ 44 

1.296 x 10“ 14 

-1.62 

-1.606 x 10“ 44 

-3.719 x 10“ 15 

-1.64 

-2.046 x 10“ 44 

2.433 x 10“ 15 

-1.66 

-2.473 x 10“ 44 

-3.020 x 10“ 15 

-1.68 

-2.886 x 10“ 44 

-1.021 x 10“ 15 

-1.70 

-3.285 x 10“ 44 

-9.549 x 10“ 15 

-1.72 

-3.670 x 10“ 44 

-4.883 x 10“ 15 

-1.74 

-4.042 x 10“ 44 

-4.188 x 10“ 15 

-1.76 

-4.400 x 10“ 44 

-2.602 x 10“ 15 

-1.78 

-4.744 x 10“ 44 

-2.001 x 10“ 15 

-1.80 

-5.075 x 10“ 44 

-1.014 x 10“ 15 

-1.82 

-5.392 x 10“ 44 

-6.463 x 10“ 15 

-1.84 

-5.696 x 10“ 44 

2.797 x 10“ 15 

-1.86 

-5.986 x 10“ 44 

-1.663 x 10“ 15 


6. Choice of Values for the Parameter n 


The parameter n controls the time evolution of the scalar field (<p). Time variations of oj, p 
and y also depend on n, as we know that both m and m 0 are functions of n. To estimate its 
value, one may use the fact that the scalar field (4>) is the reciprocal of the gravitational 
constant (G) [12, 26]. Hence, using equations (10) and (19), we can write, 

C = i = i© " = ^exp[-na((t/t o y - 1}] (54) 

According to some studies, the gravitational constant increases with time [27-29]. In the 
expression of G in equation (54), the scale factor (a) is an increasing function of time. For 
negative values of the parameter n, G would be an increasing function of time. Therefore, in 
the present study we have used only negative values of the parameter n. 


Using equations (20) and (54), we get, 
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Equation (55) suggests that the value of the parameter n can be estimated from the 
experimental observations of the fractional rate of change of the gravitational constant (G/G), 
at the present time, as reported in articles on time varying G [30]. 

According to a study by S. Weinberg, the largest possible value of (-) is 4 x 10 -10 per 

\GJ t =t 0 

year [32]. 




TABLE - 2 

Values of m 0 , m, y 0 , y D0 , (- ) for several values of the parameter 

yGJ t=t 0 

n which determines the time variation of the scalar field ( 0 ). 

n 

O ) 0 

m 

To 

Ydo 

(z c — 0.7) 

{0 

V y t=t 0 

(Tr _1 ) 

-1.63 

-1.514 

-0.183 

-1.796 

-2.566 

1.198X lO” 10 

-1.64 

-1.518 

-0.168 

-1.691 

-2.416 

1.205X 10“ 10 

-1.65 

-1.522 

-0.154 

-1.588 

-2.268 

1.212X 10~ 10 

-1.66 

-1.525 

-0.140 

-1.486 

-2.123 

1.220X lO” 10 

-1.67 

-1.529 

-0.127 

-1.386 

-1.979 

1.227X lO” 10 

-1.68 

-1.532 

-0.115 

-1.287 

-1.839 

1.234X lO” 10 

-1.69 

-1.535 

-0.104 

-1.190 

-1.700 

1.242X 10“ 10 

-1.70 

-1.538 

-0.093 

-1.095 

-1.564 

1.249X 10“ 10 

-1.71 

-1.540 

-0.083 

-1.001 

-1.430 

1.256X lO” 10 

-1.72 

-1.543 

-0.073 

-0.909 

-1.299 

1.264X lO” 10 

-1.73 

-1.545 

-0.064 

-0.819 

-1.170 

1.271X 10“ 10 

-1.74 

-1.547 

-0.055 

-0.730 

-1.043 

1.279X 10“ 10 

-1.75 

-1.549 

-0.046 

-0.643 

-0.918 

1.286X lO” 10 

-1.76 

-1.551 

-0.038 

-0.558 

-0.796 

1.293X 10“ 10 

-1.77 

-1.553 

-0.031 

-0.474 

-0.677 

1.301X 10- 10 

-1.78 

-1.554 

-0.024 

-0.392 

-0.559 

1.308X 10“ 10 

-1.79 

-1.555 

-0.017 

-0.311 

-0.444 

1.315X 10“ 10 

-1.80 

-1.557 

-0.011 

-0.232 

-0.332 

1.323X lO” 10 

-1.81 

-1.558 

-0.005 

-0.155 

-0.221 

1.330X lO” 10 

-1.82 

-1.559 

0.001 

-0.079 

-0.113 

1.337X 10“ 10 

-1.83 

-1.560 

0.006 

-0.005 

-0.008 

1.345X 10“ 10 

-1.84 

-1.561 

0.012 

0.067 

0.095 

1.352X 10“ 10 

-1.85 

-1.561 

0.016 

0.137 

0.196 

1.359X lO” 10 


Table-2 shows the values of n> 0 , m. y 0 , y D0 , (-) for different values of the parameter n 

v GJ t=t 0 

which governs the time variation of the scalar field ((f)), as per equation (10). It is found that, 
for the values of n over the range of —1.76 < n < —1.67, the values of the present EoS 
parameter (y 0 ) in this table are consistent with the range of values (approximately —1.7 < 
y 0 < — 0 . 6 ) obtained from astrophysical observations, as reported in some recent articles on 
the time variation of EoS parameter, based on an anisotropic space-time [17, 18].Over this 
entire range, the values of the parameter m are negative, indicating that o) would be less 
negative as 0 increases, as evident from equation (11), where m 0 has a negative value. For 
n < —1.81 in this table we have m > 0 , implying that oo would be more and more negative 
as 0 increases. Both of these natures of variations of <n(0) are reflected in the plots of some 
other studies on Brans-Dicke theory [13, 31]. But the second type of variation is not 
physically valid, since it corresponds to those values of EoS parameter which are not in 

conformity with the astrophysical observations in this regard. The values of (-) in this 

table are absolutely within the range of permissible values envisaged by S. Weinberg [32]. 
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Figure 1: Plots of Brans-Dicke parameter 
versus time for three different values of the 
parameter n. 



Figure 3: Plots of the equation of state 
parameter versus time for three different 
values of the parameter n. 



Figure 2: Plots of energy density versus 
time for three different values of the 
parameter n. 



Figure 4: Plot of the equation of state 
parameter versus time, for n — —1.74, on a 
larger time scale. 

















































FIGURES 
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Figure 5: Plots of density parameters, for 
matter and dark energy, for z c = 0.7, 
against the redshift parameter. 



Figure 6: Plots of density parameters, for 
matter and dark energy, as functions of 
time, for three different values of z c . 




Figure 7: Plots of the equation of state 
(EoS) parameter for dark energy versus 
time for three different values of the 
parameter n. 


Figure 8: Plots of EoS parameters versus 
time, for total energy (y) and dark energy 
( Yd )• For n = —1.74, y 0 — —0.73 and 
Ydo — ~ E04. 























































FIGURES 



Figure 9: Plots of effective interaction term 
Qeff versus time for three different values 
of z c . 



Figure 10: Plots of effective interaction 
term Q e ff versus time for three different 
values of the parameter n. 



Figure 11: Plots of effective interaction 
term Q e // versus redshift (z) for three 
different values of z c . 



Figure 12: Plots of effective interaction 
term Q e /f versus redshift (z) for three 
different values of the parameter n. 


















































7. 


Results and Discussion 


Figure 1 shows the variation of the Brans-Dicke parameter (co) as a function of time, for three 
different values of the parameter n. Here, co has a negative value and it becomes more 
negative with time. For more negative values of n, it decreases less rapidly. This nature of co, 
which is gradually becoming more negative with time, is also obtained from other studies on 
Brans-Dicke theory [13, 31]. The time dependence of co, obtained from the present study, is 
based on a scale factor which leads to a time dependent deceleration parameter (that shows a 
signature flip with time from positive to negative), unlike some recent studies [7, 14, 15] that 
are based on time independent deceleration parameters. The value of co 0 is found to be close 
to —1.55 which is consistent with the range provided in the conclusion of a recently 
published article based on Kantowski-Sachs space-time [15]. 

Figure 2 shows the variation of energy density ( p ) with time, for three different values of the 
parameter n. Here, p is found to rise steeply to a peak and then decreases at a slower rate. 
More negative values of n cause quicker attainment of the peak and also a larger peak value 
of energy density. It falls more steeply for more negative values of n. 

Figure 3 shows the time dependence of the equation of state (EoS) parameter (y) for three 
different values of the parameter n. These curves show a steep rise initially and then they 
gradually become asymptotic to a small negative value. Saturation is reached faster for less 
negative values of n. The value of y at the present time (t = t 0 ) is —0.73 which is in 
agreement with the ranges of values obtained from astrophysical observations, as reported in 
earlier studies [17, 18]. These curves are similar to those obtained from dark energy models 
based on LRS Bianchi type-V metric in the framework of Einstein’s general theory of 
relativity [18]. 

Figure 4 shows a single plot EoS parameter, for n — —1.74, on a longer time scale. It shows 
a pattern that is very much similar to the one depicted in one of the earlier studies based on 
dark energy models in anisotropic Bianchi type-I (B-I) space-time, in the framework of 
Einstein’s general theory of relativity [17]. 

The plots in the figures 2-4 show very clearly that the time evolution of energy density and 
EoS parameter, obtained here for FRW space-time in the framework of Brans-Dicke theory, 
are quite in agreement with the results obtained from anisotropic space-time in the framework 
of general relativity. 

Figure 5 shows the variations of density parameters, for dark energy and matter, as functions 
of the redshift parameter. Here, the value of p has been taken to be —0.557, causing z c — 
0.7, implying that the dark energy took over at around z = 0.7. These plots are almost 
identical in nature to the plots obtained from some recent studies [19, 20], which were carried 
out on premises completely different from the present study. 

Figure 6 shows the plots of density parameters, for dark energy and matter, as functions of 
time. Here we have three pairs (Q. m and fi D ) of density parameters, corresponding to three 
different values of the parameter p, connected to three different values of z c . 



Figure 7 shows the time variation of the EoS parameter for dark energy ( y D ) for three 
different values of the parameter n. These plots are very much similar to the ones obtained 
from studies based on anisotropic space-time in the framework of general relativity [17, 18]. 
Plots with less negative values of n take longer time to attain saturation. 

Figure 8 shows the time variations of EoS parameters for total energy and dark energy. Here 
y D has a more negative value compared with y. Their difference decreases with time, as the 
proportion of dark energy increases in the universe. For n — —1.74 the values of y 0 and y D0 
are respectively —0.73 and —1.04 respectively. These values are quite consistent with the 
ranges of values obtained from several astrophysical observations and described in some 
studies regarding dark energy and time varying EoS parameter [17-20]. 

The parameter n, in this model, plays an important role in governing the time evolution of 
various cosmological quantities. Both m 0 and m, which controls the value of co, are 
dependent upon n (eqns. 15, 16). Larger negative value of n means faster change of <p 

(= 1/G ) with time. Taking n — —1.74 we get, = 1.28 x 10 -10 Tr _1 from equation 

K - G -'t=t 0 

(55). This is well within its upper limit, i.e. 4 x 10 -10 Tr -1 , predicted by S. Weinberg [32]. 

Figure 9 shows the variation of the effective interaction term Qe// as a function of time for 
three different values of z c . It is found to decrease with time from an initial positive value, 
gradually becoming negative at a later time. This behaviour of signature flip is consistent 
with the results of a recent study on holographic dark energy [21]. For smaller values of z c 
(which correspond to more negative values of p). Q e // attains its negative phase faster. Its 
positive phase indicates that dark energy initially decayed into matter and, in the later stage, 
where Q e // < 0, matter decays into dark energy. 

Figure 10 shows the variation of Q e // as a function of time for three different values of the 
parameter n. Two of these curves decrease with time and change their signs from positive to 
negative. These two curves show that the transfer of energy, between matter and dark energy, 
reverses its direction, producing dark energy at the expense of matter at the later stage. The 
third curve, with n — —1.8, is entirely in the negative domain, becoming gradually closer to 
zero with time. It clearly shows that, there exists a value of the parameter n, below which we 
have an unidirectional energy transfer, from the matter field to the field of dark energy, at a 
gradually decreasing rate. 

Figure 11 depicts the change of Q e ff as a function of z, for three different values of z c . It is 
found that, for smaller values of z c , Q e // attains its negative phase faster. 

Figure 12 shows the variation of Q e ff as a function of the redshift parameter (z) for three 
different values of the parameter n. Like figure 10, one of these three plots also indicate the 
existence of a threshold value for n, below which the direction of energy transfer is always 
from the matter field to the field of dark energy with Q e fj becoming less negative with time. 
Above this threshold of n, the interaction term shows a signature flip from positive to 
negative as time progresses. The curve for n = —1.8 in this figure is similar to the behaviour 
obtained from a study on the interaction term between dark matter and dark energy [22]. 
Since around 80% of the matter content is dark matter, the interactions of dark energy with 
matter consist mainly of the interactions of dark energy with dark matter. 



8. Concluding Remarks 


The present study is based on a spatially flat, homogeneous and isotropic space-time, in the 
framework of Brans-Dicke theory. Here, the time evolution of Brans-Dicke parameter (m), 
energy density (p) and EoS parameter ( y ) have been determined from the field equations. 
Unlike many other studies, no cosmological constant and self interaction potential has been 
used in the formulation of the theory. The choice of a scale factor, that ensures a signature 
flip of the deceleration parameter with time, is likely to give more credibility to the results of 
this model. Without involving any self interaction potential (lA,) and cosmological constant 
(A), the time evolutions of density parameters for matter and dark energy have been 
determined and these results are sufficiently consistent with those obtained through much 
more rigorous calculations [19, 20]. The parameter p, in this model, governs the time 
dependence of Q. D , fl m and y D . It has been shown that p can be determined if one gets 
information from astrophysical observations regarding the time or redshift (z) at which we 
had = fl m in the recent past of the expanding universe. Assuming an interaction to be 
taking place between matter and dark energy, an effective interaction term has been 
calculated on the basis of the energy conservation equation. Its negative value at the present 
time implies a decay of matter into dark energy which is responsible for the accelerated 
expansion of the universe. The behaviour of this interaction term is consistent with the results 
of some recent studies [21, 22] carried out in a completely different way. The parameter n, 
which actually controls the rate of change of gravitational constant with time (as per eqn. 54), 
determines the nature of time variation of some cosmological quantities (co, p and y) 
mentioned above. It implies that the dependence of the scalar field (0 = 1/G) upon the scale 
factor (a) (eqn. 10) plays a very important role in cosmic expansion. The dependence of co 
upon cp is controlled by co 0 and m which are functions of n. Positive values of n leads to large 
negative values of y D0 , contrary to its range of values close to —1, obtained from recent 
astrophysical observations [19]. This result is in favour of using negative values of n, 
implying an increase of gravitational constant with time, as per equation (54). This behaviour 
of gravitational constant is also obtained from other studies [25-27]. Keeping in view all 
these facts, one may think of improving this model by choosing a relation, between 0 and a, 
which is different from the ansatz of equation (10), in order to have more than one tunable 
parameter like n controlling the time variation of the scalar field (0). As a future project, one 
may think of determining the time dependence of p m and p D from equations (41) and (42) by 
assuming Q to be a function of the scalar field (0) where 0 may be assumed to have a power- 
law dependence upon time, as in some earlier studies [7]. For this purpose H can be taken 
from equation (20) and y D can be assumed to have a form like y D — A yt as in equation 
(26). Here y can be replaced by y 0 , considering very slow time variation of y. as evident from 
figure 4, over a span of time close to t = t 0 . Using the solutions of equations (41) and (42), 
one can calculate Cl D and D. m , taking the total energy density (p) from equation (14). Thus it 
would be possible to establish a relationship between the density parameters (il D and fl m ) 
and the scalar field (0). 
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